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ABSTRACT 

Cold or Warm, the Dark Matter substructure spectrum must extend to objects with masses as 
low as 10 ^Mq, according to the most recent Lyman-a measurements. Around a Milky Way¬ 
like galaxy, more than a thousand of these subhaloes will not be able to form stars but are 
dense enough to survive even deep down in the potential well of their host. There, within the 
stellar halo, these dark pellets will bombard tidal streams as they travel around the Galaxy, 
causing small but recognizable damage to the stream density distribution. The detection and 
characterization of these stream ruptures will allow us to constrain the details of the subhalo- 
stream interaction. In this work, for the first time, we will demonstrate how the properties of a 
subhalo, most importantly its mass and size, can be reliably inferred from the gap it produces 
in a tidal stream. For a range of realistic observational setups, mimicking e.g. SDSS, DBS, 
Gaia and LSST data, we find that it is possible to measure the complete set of properties 
(including the phase-space coordinates during the flyby) of dark perturbers with M > 10^Mq, 
up to a Id degeneracy between the mass and velocity. 

Key words: cosmology: theory - dark matter - galaxies: haloes - galaxies: kinematics and 
dynamics - galaxies: structure - astrometry, surveys, proper motions 


1 INTRODUCTION 


Embarrassingly few observational techniques exist with which to 
reliably gauge the presence of low-mass Dark Matter (DM) clumps 
in the Universe around us. Even though swarms of completely dark 
subhaloes are predicted to orbit galaxies like our own (|Diemand| 
|et al.|[200^ [Springel et al.|[^08[ ), there is still no proven tactic 
to detect them. Until recently, the only acclaimed direct probe has 
been strong gravitational lensing (see e.g.|Dalal & Kochanek|2002| 
[Hezaveh et al.|20T4l l. However, such measurements are statistical 
in nature: only the total mass fraction locked in small sub-structure 
can be typically estimated. The largest of the small subhaloes, those 
with masses between IO^Mq and 10 ®Mq can perhaps be identified 
individually (see e.g. |Vegetti et al.|20T^|Hezaveh et al.|2013^ . Yet, 
there is little hope to prove that they do not possess a substantial 
baryonic component. 

In the Milky Way, arguably, the frontier of low-luminosity 
low-mass substructure detection was breached long ago. A popu¬ 
lation of DM-dominated dwarf satellites is known to surround the 
Galaxy (see e.g. |Aaronson|1983||Mateo|1998||Tolstoy, Hill & Tosi| 
|2009|l. For every solar mass of hydrogen, these dwarfs contain hun¬ 
dreds of solar masses of Dark Matter (see e.g. [Walker et al.|2009| >. 
Granted, the total masses of the subhaloes the dwarf spheroidals 
inhibit are not currently known. The initial mass of the faintest of 
them can be as high as 10 ®Mq (see e.g. Wheeler et al. 2015 i. How¬ 
ever, it is not implausible that they could have started with as little 


as 10 ^Mq of dm I Bland-Hawthom, Sutherland & Webster|2015| l. 
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By redshift zero, their dark masses have probably been depleted 
further while their luminosities are at the level barely seeable at 
distances far smaller than cosmological (see e.g. |Belokurov|20f^ . 

Between IO’Mq and 10 ®Mq, there ought to be roughly 2000 
DM subhaloes in a IO'^Mq host if DM were cold 
|2008| l. As of today, the census of the Milky Way’s 
tion is incomplete, but the total is predicted not to exceed several 
hundred satellites (see [Tollerud et aL]|2008[|Koposov et al.|2008| 
|2009| >. It is clear that more than half of the substructure in this 
mass regime must remain invisible. Recently, Warm Dark Matter 
has been invoked as a silver bullet to make the problems associated 
with low-mass substructure go away (see e.g. [Lovell et al.|2012| 
[2014[ l. However, in line with the up-to-date Ly a measurements, if 
the DM were warm, its mass ought to be greater than 3.3 keV at 
the 2(T level (see e.g. [Viel et al.|2013^ . This, in turn, leaves very 
little wriggle room to modify the substructure mass function above 
Vmax ~ 6 km s“' (see|Schneider et al.|2014t. 


Springel et al. 
dwarf popula- 


The next goal, therefore, is to discover subhaloes totally de¬ 
void of stars. The most promising avenue is the detection and 
characterization of the perturbations DM clumps induce on stel¬ 
lar streams in the Galaxy (e.g. |Ibata et al.|2002[[Johnston, Spergel[ 
[& Haydn|2002^ . Through a series of numerical experiments, it has 
now been established that dark subhaloes should interact with stel¬ 
lar streams sufficiently frequently, leaving small but observable sig¬ 
natures in their wake (e.g. Siegal-Gaskins&^lluriJ2008UCarlbe^ 


2009 2012[[Yoon, Johnston & Hogg[2011 Ngan & Carlberg 2014 


Ngan et al.|2015[[Carlberg|2015^ . For example, during the lifetime 
of a Pal 5-like stream roughly 5 direct impacts with subhaloes with 
IO’Mq < M < 10** Mq are expected and 20 impacts with subhaloes 
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with 10® Mq < M < 10^ Mq I Yoon, Johnston & Hogg 2011 1. There 
are several tens of tidal streams already discovered in the Milky 
Way ( [Newberg & Carlin|2016^ - plenty of potential carriers of signs 
of past interactions. 

Following the treatment introduced in |Carlberg] ( |2013) , we 
laid down the basic equations governing the formation of stream 
perturbations in |Erkal & Belokurov | p015| Paper 1). Paper I 
presents a simple picture of the three phases of stream gap evo¬ 
lution, and, most importantly, demonstrates that a measurement of 
the gap density profile alone is not sufficient to break the strong de¬ 
generacy between the perturber’s mass and size, the impact geom¬ 
etry and kinematics, and the time since flyby. However, if comple¬ 
mentary phase-space measurements around the gap were collected, 
would the evidence assembled be ample to reconstruct what ex¬ 
actly happened during the run-in? This is the question we strive to 
answer in this Paper. 

In this work we will simulate cold streams similar to that of 
Pal 5 (see e.g. |Odenkirchen et al.|2003| l. We will consider the 6d 
phase-space structure of the gap and show that using this expanded 
set of observables, it is possible to infer the properties of the sub¬ 
halo’s flyby, up to a Id degeneracy. This is a degeneracy between 
the subhalo mass and the velocity of its flyby which is unavoidable 
since it is actually a degeneracy in the velocity kick imparted on 
the stream particles during the flyby. We will show that with proper 
motions and distances from Gaia, Dark Energy Survey (DES)-like 
photometry, and moderate resolution spectroscopy for radial veloc¬ 
ities, it will be possible to detect and characterize dark subhaloes 
down to M ~ 10 ^Mq within 10 kpc. We will also discuss the mass 
range which can be probed by upcoming surveys like the Large 
Synoptic Survey Telescope (LSST) and give general expressions 
for the mass sensitivity as a function of the errors in each of the 6d 
phase-space coordinates. 

Finally, we deliberate on the constraints that can be derived if 
only a subset of the observations can be secured. For example, for 
sufficiently distant streams it may not be possible to get proper mo¬ 
tions with Gaia. We show that in general one only needs to measure 
the following: the shape of the gap on the sky, the density profile of 
the gap, and the line-of-sight velocity along the gap. Fortunately, 
these are the simplest measurements to make for targets across a 
large range of heliocentric distances, so this does not limit our abil¬ 
ity to infer subhalo properties. 

This paper is organized as follows. We give the setup and 
motivation in Section]^ In Section]^ we present analytic expres¬ 
sions for the 6d structure of the gap. Next, in Section]^ we use 
these results to infer the properties of a subhalo from a gap in 
an N-body simulation without accounting for observational errors. 
In Section]^ we discuss the information content of each observ¬ 
able and find the minimal set of observables which are needed to 
constrain subhalo properties. In Section]^ we infer the properties 
of subhaloes from gaps in simulations using observational errors 
of current and future surveys. Continuing with realistic errors, we 
present analytic expressions for the mass range accessible as a func¬ 
tion of survey errors in Section]^ In Sectionj^we discuss complica¬ 
tions which will arise when extending this model to more realistic 
streams. Finally, we conclude in Section|^ 


2 SETUP AND MOTIVATION 

The aim of this work is to show that it is possible to start with a 
gap in a tidal stream and infer the properties of the subhalo which 
created it. In order to study this problem, we will consider a simple 


setup with a stream on a circular orbit which is perturbed by the 
flyby of a Plummer sphere. This is identical to the setup of Paper I. 
In that work, we were primarily interested in the density structure 
of the gap and we showed that a measurement of the stream density 
profile near the gap does not uniquely constrain the properties of the 
subhalo which created the gap. In fact, it only provides constraints 
on three out of the seven parameters which describe the subhalo 
flyby. Thus, more constraints are needed. The results of Paper I 
also contain predictions for the entire 6d structure of the gap which 
we will exploit in this work. With this extra information, we will 
show that we can constrain the subhalo properties down to a Id 
degeneracy between mass and velocity. 

One of the main benefits of this simplistic model for gaps is 
that all of the predictions are analytic. With these expressions for 
the stream shape, we can predict the subsequent evolution of the 
stream after a subhalo flyby. For a given gap, the model can then be 
used to infer the properties of the subhalo which created it. Further¬ 
more, the analytic expressions will allow us to probe the constrain¬ 
ing power of each observable and we will discuss what degenera¬ 
cies remain if only certain observations are possible. Finally, we 
will use these expressions to give constraints on what mass range 
is gaugeable given the errors of current and future surveys. Without 
these analytic predictions, a lot of this would not be possible and we 
would need to build intuition entirely through individual numerical 
examples. 

Figure [T] displays the geometry of the subhalo-stream en¬ 
counter as used in this work. The stream is on a circular orbit and 
we choose our coordinates to be centered on the point of closest 
approach along the stream with x denoting the radial direction, y 
denoting the direction along the stream, and z denoting the direc¬ 
tion perpendicular to the stream plane. This is identical to the setup 
in Paper I. We assume that the stream is orbiting in a spherical po¬ 
tential which is otherwise unconstrained. The stream is at a radius r 
with a circular velocity Vy = yjrdr4>- This model allows for any im¬ 
pact geometry. We denote the subhalo velocity as {wx,Wy,w,) and 
the stream velocity as (0, Vy, 0). At the point of closest approach, the 
vector between the stream and the subhalo must be perpendicular 
to the relative velocity {Wx, Wy — Vy,w,). This relative velocity, along 
with the impact parameter b, defines the subhalo’s trajectory up to 
an overall sign since the subhalo could pass on either side of the 
stream. We will encode this sign in the impact parameter, b. Next 
we have the mass of the subhalo, M, and its scale radius, r^. Fi¬ 
nally, we have the time since impact, t. For a given gap, we take the 
impact point to be the current position of the gap center integrated 
backwards in time by t. Thus we have seven model parameters: 
three velocities, v?, the impact parameter, b, the mass M, the scale 
radius, r,, and the time since impact, t. Following the convention 
of Paper I, we break the relative velocity between the stream and 
subhalo into a velocity along the stream, W|| = Vy—Wy, and a perpen¬ 
dicular velocity of w± = ■\/w^ + w^. The magnitude of the relative 

velocity is denoted by w = -l- w^. 

During the flyby, we assume that the stream is moving along a 
straight line. This is needed to compute the velocity kicks along the 
stream. Since the velocity kicks take place over a relatively small 
distance, several kpc at most for the most massive subhaloes we 
will consider, this is a reasonable assumption. 


3 NEW OBSERVABLES 

The first observable we will consider is the density profile along the 
stream as in Paper 1. In addition, we will now use the analytic model 
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Figure 1. Our axis convention looking down on the stream. The dotted line 
is the stream and the dashed line is the path of the suhhalo. The solid lines 
show the impact parameter, b, and its angle in the x-z plane, a. The center 
of the host potential is to the left and the orbit is counterclockwise. Note that 
the angle a can he re-expressed in terms of the velocities Wj^ and since the 
velocity is perpendicular to the impact parameter at closest approach. This 
determines the impact parameter up to an overall sign which we include in 
b. We note that this figure is identical to Figure 1 in Paper I. 


produced in Paper I to make predictions for the 6d phase space 
structure of the stream after the flyby, giving five new observables. 
First we have the spatial shape of the stream near the gap which 
gives two observables: 1) the distance to the stream as a function of 
the angle on the sky, xwsil/ (where xj/ is the angle along the stream) 
and 2) the shape of the stream on the sky, z vs if/. Next we have the 
velocity structure of the gap which gives three observables: 3) the 
radial velocity as a function of the angle on the sky, x vs if/, 4) the 
proper motion along the stream as a function of the angle on the 
sky, y vs if/, and 5) the velocity perpendicular to the orbital plane 
as a function of the angle on the sky, z vs if/. We note that these 
observables are what an observer in the center of the galaxy would 
measure. They do not account for the solar position or velocity. 
We will give analytic predictions of these observables below and 
then use them to infer the subhalo from several N-body simulations. 
Note that these predictions can only be expressed parametrically, 
i.e. to describe the shape of the gap on the sky (z vs if/) at time t, we 
express both z and if/ in terms of the angle from the point of closest 
approach, if/Q. 

Crucially, the structure of the gap is controlled by the velocity 
kicks imparted on the stream during the flyby. The kicks in the 
x, y, and z direction (radial, along the stream, and perpendicular 
to the orbital plane), as a function of distance from the point of 
closest approach were computed analytically for a Plummer sphere 
perturber in Paper I: 


AVx(iAo) = 


2GM ^ - ^I'oroWiiW,) 


t^Vy{lf/0) = ■ 


'IGMxj/Q 


0 


) 




wro 


2GM [bw^^+doroW\\w^ 
Av.((Ao) = - 


{K- 




( 1 ) 


These formulae are slightly different than the ones in Paper 1 since 
we have re-expressed the kicks in terms of the three components 
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of the flyby velocity instead of two velocity components and an 
angle. We have also re-expressed the distance from the point of 
closest approach, y, in terms of the angle along the stream, i.e. y = 
XQif/Q. Next, we list the prediction for each of the 6d phase space 
coordinates. 


3.1 if/ 

First we have the angle along the stream. This is given in equation 
16 of Paper I: 


iA(iAo. 0 - </'o + 


f<ko-g 
if/l + B^- 


( 2 ) 


where if/o = denotes the initial distance between a point along 
the stream and the point of closest approach, 

4-ft 4sin(^) 2(1 - cos(-;^)) n,||Vr;t ro 


f = 


7 


7 


7 / 


(3) 


2{l-cosC-^)) ro 

7^ r(,w\ VyP 


(4) 


B- = 


b^- 


(5) 



and the timescale t is given by 


T = 


2GM 


( 6 ) 


(7) 


Note that these expressions are slightly different from those in Pa¬ 
per I since we have re-written the angle a in terms of ratios of Wj,, 
Wj., and w±. The density is related to this angle sincep oc I 


3.2 x 

The radial distance of each particle in the stream is computed in 
equation 8 of Paper I: 

^ 2roAv,(<Ao)(l-“s(^)) ,„Av,(^^o)sin(^) 

Ax{^o,t) =-5-+-, (8) 

Vy 7^ Vy y 

where the velocity kicks, AV;,.(if/o), Avy(if/o), are specified in 0. 


3.3 z 

The kicks in the z direction tilt the orbital plane of the resulting 
orbits by an angle of Av^/v^, in the x direction, therefore the offset 
in stream latitude of the particles is given by 


AS^(if/o, t) 


Av^fiAo) 


Vyt 

sm( —), 
''0 


(9) 


where AS, = — which denotes the difference in stream latitude 

ro 

from the unperturbed stream track. 
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3.4 i: 

The radial velocity can be derived from a time derivative of § 
which gives 


Ax(i/fo, t) 


2 Avv(i/'o) yvvf yvJ 

-sm( —) + Av^(^/^o)cos( —). 

7 ra ro 


( 10 ) 


3.5 y 

The velocity along the stream is given by 


Ay(iAo, t) = (ro + Ax(iAo, f)) e(if/o, t)-Vy, (11) 


where ^(i/'o, t) is the angular velocity along the stream and is given 
in equation 12 of Paper I and Ax(tl/Q, t) is given in Plugging in 
these expressions and taking only the leading order terms we get 


Ay(if/o,t) = -Avyiif/o)- 


■ 7 ^ cos(2^) 

^+2Av,(<Ao)- 


sin(^) 

-AvA^ka) -^ 

7 

( 12 ) 


We note that we only consider the leading order terms in Avj: and 
Avj, in this equation since the analysis in Paper I was done at leading 
order. 


3.6 z 

The velocity in the z direction comes from a time derivative of the 
offset in the z direction, i.e. the time derivative of 

Az(iAo, t) = AvA^Ja) cos( —). (13) 

ro 

By combining these six results, and we have a para¬ 

metric model for the phase space structure of the stream at time t 
in terms of the initial position along the stream, if/Q. In Figure]^ 
Figure]^ we show a comparison of the model presented in this 
section against the results of N-body simulations which we will de¬ 
scribe in the next section. From these Figures, it is evident that the 
model can successfully reproduce the structure of the gap. 


4 INFERENCE WITHOUT ERRORS 

Now that we have introduced the new observables, we will use a 
suite of N-body simulations to demonstrate that the model works. 
We will then run a set of MCMCs to show that the model can be 
used to constrain the properties of the subhalo which created the 

gap- 


4.1 Simulations 

To generate the gaps we carry out three N-body simulations with 
the pure N-body part of gadget-3 which is closely related to gadget- 
2 ( |Springel|2005^ . 

The setup is similar to the realistic simulations in Paper I. We 
first generate a stream by taking a globular cluster modelled as a 
King profile with a mass of 10 ^Mq, a core radius of 13 pc, and 
w = 5. This progenitor is placed on a circular orbit with a radius of 
10 kpc around an NFW potential ( [Navarro, Frenk & White|19^ 
with M = lO'^Mo, c = 15 and = 14.0 kpc. The circular velocity 
at 10 kpc is Vy = 168.2 km/s. The King profile is represented with 
10® particles which have a smoothing length of 1 pc. The cluster is 


evolved for 5 Gyr and in this time a cold stream with a total length 
of ~ 180° is produced. We take this stream and place a subhalo on 
an orbit to impact the center of the leading arm of the stream, which 
corresponds to a radius of 9.80 kpc, at 150 km/s perpendicular to 
the orbital plane. We then follow the evolution of the stream for an 
additional 5 Gyr. We have chosen a simple impact geometry in this 
example but the analytic model works for any impact geometry. 

In our experiments, we use three different subhaloes to show¬ 
case the range of masses and scale radii for which we will be able to 
infer the subhalo properties. These are chosen to roughly lie in the 
center of the scale radius-mass relation in the public data release of 
the Via Lactea II simulation jPiemand et al.|2008t . First, we have 
the largest subhalo with M = 10*Mq with = 625 pc. Next, we 
have an intermediate subhalo with M = IO^^Mq with r, = 395.3 
pc. Finally, we have the smallest subhalo with M = irPM q and 
r, = 250 pc. We note that while these subhaloes were chosen in 
the middle of the scale radius-mass relation from Via Lactea II, we 
have modelled them with Plummer spheres instead of NFW pro¬ 
files. As a result, our Plummer spheres reach a maximum circular 
velocity at a smaller radius than if they were modelled with NFW 
profiles. If we instead compare the radius at which the maximum 
circular velocity of our subhaloes is reached to those in Via Lactea 
II, we find that our subhaloes are still within the scatter but are now 
closer to the edge of the relation rather than the center. 

As we saw in Section]^ the amplitude of the oscillations in¬ 
duced by the subhalo’s passage is governed by the maximum veloc¬ 
ity kick induced by the subhalo. For a direct impact, this maximum 
is proportional to M/r^ which follows from 0. The properties of 
these subhaloes are arranged so that the size of the perturbation 
decreases by a factor of two as we go from the most massive to 
the intermediate subhalo, and another factor of two as we go from 
the intermediate to the least massive subhalo. Finally, we have to 
choose how long after impact to look at each of the gaps. We choose 
to look at the perturbations when the deflection on the sky (z vs tj/) 
is relatively large at t = 450 Myr. Note that the different subhaloes 
will have created different sized gaps by this time. While the gap 
will grow in time, becoming easier to observe, we choose to look 
at the gap after a short time where our model is most accurate. At 
late times, the distribution of energy and angular momenta along 
the stream can cause the gap to fill in which is not captured by our 
model. We choose to look at all subhaloes at the same time to make 
it a fair comparison. 

In Figure [^3|4| we show the six observables for each of the 
three instances of the perturbed stream at these times. Note that 
the scales on the axes are different in each of the figures. The grey 
bands in these figures show the width of the stream and the black 
curves show the prediction from the model for each setup. We stress 
that the black curves are not fits to the data but rather a prediction of 
the model using the true parameters. While the model agrees well 
with the simulations, it does not capture the entire behavior. For 
example, if we look at the v,,. panel in Figure]^ or Figure]^ we see 
that the model does not pass through the center of the grey band. 
This slight disagreement is important as it will result in biases when 
we fit our model to the simulations as demonstrated in Section [4] 
and Sectionj^ 

In these simulations, we have focused on direct impacts which 
will produce the largest signatures. The analytic results of Section|^ 
allow us to explore how large the signatures will be if we change 
the impact parameter, the mass of the subhalo, or the velocity of 
the flyby. The size of the perturbation is governed by the size of 
the velocity kick which can also be estimated qualitatively by tak¬ 
ing the typical acceleration during the flyby, multiplied 
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Figure 2. Gap observables for a flyby of a Plummer sphere with a mass 
of lO^M© a scale radius of = 625 pc at r = 450 Myr. The grey bands 
shows the Icr scatter of the N-body simulations and the black curve shows 
the curve predicted by the analytic model. Note that the scales on the y axes 
are different from Fig.[^and Fig.|^ 



V'(°) 


Figure 4. Gap observables for a flyby of a Plummer sphere with a mass 
of lO^M© a scale radius of = 250 pc at r = 450 Myr. The grey bands 
shows the Icr scatter of the N-body simulations and the black curve shows 
the curve predicted by the analytic model. Note that the scales on the y axes 
are different from Fig. Hand Fig.[^ 



Figure 3. Gap obsei'vables for a flyby of a Plummer sphere with a mass 
of 10^ a scale radius of = 395 pc at t = 450 Myr. The grey bands 
shows the Icr scatter of the N-body simulations and the black curve shows 
the curve predicted by the analytic model. Note that the scales on the y axes 
are different from Fig.[^and Fig.|^ 


■i/i,2+r2 

by the duration of the flyby, ~ ' , ‘ , to give a velocity kick of 

- ^ . As we showed in Figure 3 of Paper I, the maximum kick 

in the y direction is given by so the qualitative approach 

almost gives the correct behavior. Thus we see that as we increase 
the mass of the subhalo, we will get a larger signature. Likewise, if 
we increase the velocity of the flyby, the impact parameter, or the 
subhalo’s size, we will get a smaller signature. The maximum size 
of the velocity kicks in the r and z can also be computed using 0 
and exhibit roughly the same behavior although the dependence on 
the velocities, r^, and b is more complicated. 


4.2 Inference 


We will now infer the properties of the subhalo without accounting 
for any observational error. This showcases how well this method 
works in the best case scenario. In Section we will do the same 
exercise accounting for observational errors for both contemporary 
and future surveys. We use emcee jForeman-Mackey et al.|20T^ 
to explore the likelihood landscape of our models. 

The MCMC samples the seven model parameters, 
M,r,,Wx,Wy,w^,b,t, as well as 6 nuisance parameters for 
overall shifts in the 6 phase space coordinates. These nuisance 
parameters are needed to account for possible offsets of the center 
of the gap due to epicyclic motion. As we can see from l[TJ, the 
center of the gap receives no kick along the stream so it has the 
same orbital period as the unperturbed stream. However, it can 
receive kicks in the other directions which causes the center of 
the gap to execute epicycles around its unperturbed position. As a 
result, when we measure the location of the gap, we do not know 
the current offset of the center of the gap from the unperturbed 
stream and we need to marginalize over this. 

The model described in Section gives the centroid of the 
gap in all six observables. The centroid of the gap in the N-body 
simulations is found by fitting a Gaussian profile (controlled by 2 
parameters describing its center and intrinsic width) to each ob¬ 
servable in 1 degree bins along the stream. This gives the centroid 
as well as its error for all observables except the density which is 
given Poisson errors. This procedure mimics the process of decou¬ 
pling the stream and the background, where, at least for some of 
the observables, the separation is done in bins of angle along the 
stream and not per star. 

With regards to modelling the subhalo-stream encounter itself, 
the likelihood for each observable is simply: 



ln(2;rcr^) 

2 


(o, - m,-)^ \ 

2<t^ I’ 


(14) 


where i denotes the angle bins along the stream, o,, cr, are the cen¬ 
troid and its error in the N-body simulation, and m, is the centroid 
from the model. We assume the observables are independent so the 
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full log-likelihood is the sum of the individual log-likelihoods for 
each observable. 

In the analysis presented here, our priors are as follows. For 
the mass, we assume a uniform prior in logM between a mass 
of 10 ^Mq and 10 ®Mq. For the scale radius, we assume a uniform 
prior between 1 pc and 2 kpc. For the impact parameter, we take a 
uniform prior between -2 kpc and 2 kpc. The prior on the veloc¬ 
ity of the subhalo is a Maxwellian distribution with a dispersion, 
cr, related to the circular velocity at 10 kpc, Vcirc, by the relation 
O' - Vcirc/ Vs. For our setup this gives cr = 97.1 km/s. This prior 
is very important since, as we will see and discuss below, there is a 
degeneracy between the mass and velocity of the flyby. For the time 
since impact, we assume a uniform prior between 100 Myr and 2 
Gyr. Next we have priors for the nuisance parameters. The shift in 
the angle along the stream is uniform from -5 to 5°. The shift in 
the radial distance is uniform from -0.2 kpc to 0.2 kpc. The shift 
in all velocity components is uniform from -2 to 2 km/s. 

In Figure]^ we show the result of the MCMC modelling for 
the lowest mass subhalo. We see that without accounting for ob¬ 
servational errors, this method can accurately infer the properties 
of subhalo with a mass of IO^Mq. Almost all of the parameters are 
well constrained. The only exception is the degeneracy between the 
mass and the velocity which we will discuss in detail in Section]^ 
This degeneracy is due to the fact that the model cannot infer M and 
w independently, but rather can only deduce their ratio M/w. As a 
result, the constraint on the mass depends on the prior we place on 
the subhalo velocity. Note that the likelihood does not peak at the 
true value for every parameter. In fact, Wj., b, and especially t are 
discrepant. This is likely due to the fact that our model does not 
account for all of the features in the stream, such as epicyclic feath¬ 
ering (e.g. |Kupper, MacLeod & Heggie|2008[|Kupper et al.|2010^ 
and that we assume that the stream is on a circular orbit. While the 
progenitor itself is on a circular orbit, the stream it generates will 
be slightly eccentric. For the subhaloes with and 10 *Mq, 

the figures look similar albeit with slightly less scatter. 


which are a function of ijjf,. If any of these coefficients are differ¬ 
ent, a different gap will be produced. In the ideal case, we would 
measure these functions directly and then read off the coefficients. 
Degeneracies can only arise if there are directions we can move 
in the 6d parameter space, p, which do not change any of these 
coefficients. If we consider a small perturbation in our parameter 
space, Sp, the degeneracies would need to satisfy a system of linear 
equations: 

Sp-WAi = 0, (16) 

where the gradient is taken with respect to p. This is a system of 
6 equations with 6 unknowns. In principle, this could have one so¬ 
lution and no degeneracy. However, for generic parameters, these 
6 constraints are not linearly independent and give rise to only 5 
independent constraints. Therefore, there is a Id degeneracy which 
involves a scaling of mass and velocity: 

M ^ AM 

Wjc -» AWjc 
W|| —> Awii 

w, Aw,. (17) 

This degeneracy is insurmountable since it does not affect the kicks 
themselves and hence has no effect on the gap. 

5.2 Degeneracies in the observables 

We can now apply this same argument to the observables described 
in Section |3l We can then combine various subsets of these obser¬ 
vations to see what these combinations can constrain. The best we 
can hope for is to match the ideal case of a Id degeneracy. The ob¬ 
servables have an additional parameter which must be fit: the time 
since impact, t. Thus when considering observables we have a 7d 
parameter space which we will denote, q = [M, Wx, Wy, w„ r,, b, 1], 


5 DEGENERACY AND INEORMATION CONTENT 

As we saw in Figure]^ there is a degeneracy between the mass and 
velocity of the subhalo. We will now study the set of degeneracies 
that arise when inferring stream properties. 


5.1 Degeneracy in the velocity kicks 


First, we will study the degeneracy present in the velocity kicks, 
AVj-, AVj,, AVj, induced by the subhalo’s passage. These kicks gov¬ 
ern the subsequent evolution of the gap so the best we can do is 
recover the information content encoded in the kicks. These kicks 
are given in (jTJ, which we will re-write as 




Avyipo) = 


Av^fiAo) = 


A lip) + i/'oA2(/0 
lAo + APp) 

>kl+A.3iP)' 
APp) + ifJoAPp) 
lAo + APP) 


(15) 


where we have re-expressed in terms of the six coefficients, 
Aj, which depend only on the parameters p = {M,Wx,Wy,w,,rx,b], 
which describe the flyby. For completeness we give the A, in ([ATJ. 
The entire effect of the flyby is encapsulated in these three kicks 


Gap density 


The first observable is the gap density along the stream. This den¬ 
sity is derived from the map between initial and final conditions 
and is given in We re-write this in terms of the notation in 0: 


t//(po, t) = tf/o + 


Blip) + ipoBiiP) 

i/'o + Bpp) 


(18) 


where the coefficients Bj are only functions of p and are indepen¬ 
dent of ij/f,. 


X VJ ijjo 


Next we have the radial distance along the stream, which we 
re-write as 


AxipoO) 


Bpp) + poBsiP) 

4'1 + Biip) 


(19) 


Z Vi l//o 

The stream latitude is governed by l|^, which we re-write as 


Az(iAo,0 


B6iP) + ipoBiip) 

if/l -H BPP) 


( 20 ) 
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M(IO^Mq) (km/s) Wj, (km/s) w, (km/s) r, (kpc) b (kpc) t (Myr) 


Figure 5. Result of MCMC for 10^ Mq subhalo with = 250 pc with no errors. Almost all of the subhalo properties are well constrained. However, there is a 
degeneracy between the mass and the velocity which we will discuss in Section]^ 


X vs 


zvsxjjo 


The radial velocity is given by jl0|l, which we re-write as 




Bn(4) + 


( 21 ) 


The proper motion in the z direction is give in ( |13^ , which we re¬ 
write as 


Az((Ao, t) 


'Pi + 


(23) 


y vs if/Q 


The proper motion along the stream, 6, is given in 1 12 i, which we 
re-write as 


AyiiAo, t) 


B[o(q) + ipoBii(d) 

tpl + Biiip 


( 22 ) 


Note that the denominator is identical for all of these observables. 
This is because all of these terms are related to the velocity kicks 
imparted by the subhalo which have the same denominator (see the 
denominator of |T](). For completeness we give the B, in ( |A2^ . 

In the expressions for the observables above, jl8| > - ( |23^ , there 
are 13 separate coefficients which govern the observables. If all 
observables are measured, we can determine the resulting degen¬ 
eracies in the parameters, by trying to find variations of the pa- 
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rameters, i.e. S4, which do not change the i?, : 

64-WBi = 0, (24) 

where the gradient is taken with respect to q. By carrying out this 
procedure, we find that these 13 constraints on 7 unknowns give rise 
to 6 linearly independent constraints, leaving the same degeneracy 
as above, i.e. |T7}. This means that if all observables are measured, 
we can recover the same information as if we measured the kicks 
themselves. 

Note that this result is for a generic value of the parameters. 
For specific values of the parameters, it is possible to have addi¬ 
tional degeneracies if some of the coefficients vanish, i.e. if S, = 0 
for some i. For example, since the observables oscillate in time, 
there will be times when some of the observables are flat and con¬ 
tain less information. Another example is if = 0. In this case, 
as we can see from the velocity kicks in l|TJ, Av^, and AVj only de¬ 
pend on the combination -l- rj and not b or r, independently. As 
a result, if the observables do not constrain (or poorly constrain) 
the shape in x and x, we should expect an additional degeneracy 
to appear between b and where b- + is constant. This is not a 
true degeneracy since the gap size does not only depend on E + rl. 
However, since the leading order term in the gap size only depends 
on b^ + r], the gap size is not very sensitive to b. As a result, b and 
r, can be effectively degenerate. We will see an example of this in 
Section IbAl 


5.3 Degeneracies in typical observational setups 

Now that we have discussed how to classify the degeneracy when 
using all observables, we will consider the degeneracies present in 
common observational setups which only use certain observables. 
We will repeat the procedure above, using only the coefficients 
present in the selected observables. 


5.3.1 Gap density, gap shape in z, and radial velocity 

One of the simplest useful setups involves measuring the density of 
the stream on the sky, the shape of the stream on the sky, and the 
radial velocity profile along the stream. These three observables de¬ 
pend on the coefficients fii,B 2 .S 3 . Se. S 7 ,Ss. and 69 . Using these 
coefficients, we assemble and solve the system of 7 linear equations 
as in Interestingly, these 7 equations also result in 6 linearly 
independent constraints on the 7d parameter space, 64, leaving the 
degeneracy between mass and velocity, i.e. 0 - This setup is quite 
simple since it only requires the radial component of the velocity 
which is far easier to measure than proper motion. This approach 
should also allow us to extend the search for gaps out to larger 
distances since the errors for these observables do not depend as 
sensitively on distance as the errors in proper motion. 


5.3.2 Gap density and gap shape in z 

Another approach is to measure the gap density and gap shape on 
the sky which can be computed using only photometric surveys and 
do not require any velocity information. These observables depend 
on the coefficients Si, S 2 , S 3 , Sg, and S 7 . We setup and solve the 5 
linear equations as in \2A) and find 5 independent constraints, leav¬ 
ing a 2d degeneracy. One of the dimensions of this degeneracy is 
the mass-velocity degeneracy discussed above. The second degen¬ 
eracy is more complicated and is between all of the parameters. 


5.3.3 Gap shape, radial velocity, proper motions 

In this approach, we make use of all of the observables except the 
density along the stream, i.e. x, z, x, y, and z vs [j/. This setup makes 
use of all the coefficients except B[ and 62 - We setup and solve the 
11 linear equations as in {24) and find 6 independent constraints 
leaving the degeneracy between mass and velocity. This is a useful 
setup since we need to understand the stripping rate of the pro¬ 
genitor as a function of time to fully model the density along the 
stream. Since this may prove challenging, it is useful to know that 
the density along the stream is not needed to constrain the subhalo 
properties. 

In Section [ 6 ^ we will compare the inferred subhalo properties 
when using these three subsets of observables against the inferrence 
when using all observables. 


6 INFERENCE WITH IMPERFECT DATA 

Now we will repeat the inference done in Section [4^ and account 
for observational errors. 


6.1 Observational Errors for Current and Future Surveys 


We will now assume three different observational setups to esti¬ 
mate the subhalo mass range accessible with this technique. Their 
details are available in Table [T] The first setup we consider fol¬ 
lows closely the properties of the SDSS ( |Ahn et ^|2012[ ) data. 
Next, we have the “Super-SDSS” survey which consists of DES- 
like ( [The Dark Energy S urv ey Collaboration|[2005| l photometry, 
coupled with Gaia 1 Perryman et al.||2001^ astrometry for proper 
motions, and deep multi-object spectroscopy from VET (e.g. Ko- 
posov et al.||2011| |Gilmore et al.||2012J, WEAVE (|Dalton et al. 
2012 | >, and 4MOSTlfe Jong et al'|2012^ for radial velocities. Fi¬ 
nally, we have the LSST setup, which is similar to Super-SDSS 
but uses LSST ^LSST Science Collaboration et al.|2009^ photom¬ 
etry and Gaia-)-LSST for proper motions to extend observations to 
much fainter stars. 


6.2 Realistic Stellar Observables 

In order to understand the subhalo mass range which can be probed 
by these surveys, we need to populate the stream with a realistic 
distribution of stars. We use the N-body simulations described in 
Section [4T] which model the globular cluster with identical mass 
particles. We use PARSEC stellar evolution models ( [Bressan et al.| 
|2012[[Chen et al.|2014[|Tang et al.|2014| l to assign masses and lumi¬ 
nosities to each particle. We chose the initial mass function (IMF) 
from |Chabrier| ( ( 200 T| > and used a single isochrone with an age of 10 
Gyr, and a metallicity of 0.01 Zg. We made use of two photometric 
systems: SDSS ugriz and Gaia G. The stellar mass of each particle 
was sampled from the corresponding present day mass function un¬ 
til the total stellar mass was equal to 10 ^ Mg, the mass of the cluster 
used in the simulation. This resulted in approximately 360,000 N- 
body particles being tagged as stars. For each star, the magnitude in 
each band was then linearly interpolated from the tabulated values 
in the isochrone table. Given the magnitudes (g, r, G) of each par¬ 
ticle, observational noise was added to the position, velocity, and 
magnitudes using the errors and cuts in Table [T] We note that the 
number of stars above the r-band cutoff in Table [T] and within 20 
degrees of the gap center increases dramatically as we consider the 
deeper surveys. The SDSS setup has 1404 stars, the Super-SDSS 
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Setup 

a, 6 

RV 

Proper Motion 

Distance 

SDSS 

Super-SDSS 

LSST 

0.1”, Oim = 22 

0.1”, Him = 24 
0.1”, = 27 

5 km/s at r = 16, 20km/s at r = 20 

1 km/s at r = 16, 5 km/s at r = 21 

1 km/s at r = 16, 5 km/s at r = 21 

None 

Gaia 

Gaia+LSST 

None 

DES photometric accuracy 
LSST photometric accuracy 


Table 1. Observational uncertainties for contemporary and upcoming surveys. The first column gives the name of the observational setup. The second column, 
a, 6, gives the observational error and limiting r-band magnitude for determining positions on the sky. The third column, RV, gives the uncertainty in deter¬ 
mining the radial velocity as a function of magnitude. We linearly interpolate the errors between the bounds given. The fourth column. Proper Motion, gives 
the proper motion errors. The Gaia proper motion errors are given in the Astrometric Performance section of http://www.cosmos.esa.int/web/gaia/science- 
performance The LSST proper motion errors are given in Table 3.3 of |,LSST Science Collaboration et al!]j2009) . The fifth column. Distance, gives the source 
of distance en'ors. The distances are calculated using photometric distances as described in Sec.|6.2|and distance enors come from the photometric accuracy 
of each survey. The DES photometric accuracy is discussed in the text. For the LSST photometric accuracy we use the fits given in Table 3.2 of|LSST Science] 
ICollaboration et ^j2009) . 


setup has 4039 stars, and the LSST setup has 10735 stars. We stress 
that we do not account for contamination from Milky Way stars or 
from background galaxies. 

The relative distances along the stream were calculated us¬ 
ing the photometric parallax method. Given the isochrone of the 
stream, we use the g - r color to infer the r-band absolute magni¬ 
tude. The apparent r-band magnitude combined with the absolute 
magnitude gives the distance modulus and hence the distance. We 
propagate the errors from the photometric accuracy for each color 
and use the slope of the isochrone to give the error in absolute mag¬ 
nitude. We do not account for uncertainties in the isochrone itself 
as we are primarily interested in relative distance changes across 
the gap. For the SDSS setup, we chose not to use any distance 
measurements since the distance errors are too large to give use¬ 
ful constraints. For the Super-SDSS setup, we use the expected 
photometric accuracy of DES. To estimate the single epoch pho¬ 
tometric accuracy, we use the photometric catalog described in |Ko-| 
Iposov et^ ( |2015| l which is derived from publicly available images 
from the first year of DES. We assume that all five epochs of DES 
will have a similar accuracy to get a final photometric accuracy 
which is V5 smaller than the single epoch accuracy. Eor complete¬ 
ness, we list our expected five epoch DES r-band magnitude errors: 
0.0005, 0.0096, 0.0019, 0.0038, 0.0079, 0.017, 0.039, 0.091 for r- 
band magnitudes of 17,18,19,20,21,22,23,24 respectively. The g- 
band magnitude errors are similar to these. For the LSST setup, we 
use the LSST photometric accuracies which are reported in |LSST| 
Science Collaboration et al.| ( |2009| l. 

6.3 MCMC Details 

The details of the MCMC are almost identical to those presented in 
Section !?!^ The only difference is that the observables for each star 
are now scattered by the appropriate observational error and this is 
accounted for when we compute the centroid and its error for each 
observable along the stream. 

In Figure [6l7|8| we show the six observables and their Icr er¬ 
ror bars for the W’Mq subhalo for the three different observational 
setups: LSST, Super-SDSS, and SDSS. We note that the error bars 
in these figures are the errors on the mean within each angular bin 
and are therefore significantly smaller than the error on any individ¬ 
ual measurement. We see that, as expected, the shape of the stream 
on the sky and the density have the smallest error bars, followed 
by the distance along the stream, and then the radial velocity. We 
see that the error bars for the Super-SDSS setup are already quite 
constraining. As we move from Super-SDSS to LSST, the main 
improvement is the superior photometry and depth (i.e. number of 
stars) which in combination give smaller error bars in the shape 



Figure 6. Observables for 10 ^Mq subhalo with LSST errors. The red curves 
show the maximum likelihood result and the blue error bars show the Icr 
errors of the mock observation. 

of the stream gap on the sky as well as its density. There is also a 
noticeable improvement in the proper motion measurements. 


6.4 Inference 

Figure 1^ shows the Icr constraints on M, r^, b, and t for each sub¬ 
halo and each observational setup. For and t we see that the 
constraints get tighter as we use higher quality data. The exception 
is M where the constraints do not improve with better data. This is 
due to the degeneracy between M and w which means that the con¬ 
straint on M comes mostly from the prior on the subhalo velocity. 

Figure[^gives the Icr constraints on M and rj for the 10 */Wq 
subhalo using LSST errors. In this figure, we vary the observables 
used in the MCMC. We see that if we use all of the information, or 
subsets of the observables which can constrain all of the properties, 
we obtain results of similar quality. However, if we have too few 
constraints then the error bars get larger. The three subsets we use 
here are described in Section [53] 

In Table we give the 5-95% confidence intervals for M and 
r, for each subhalo and observational setup. The trends here are 
identical to those seen in Figure]^ In Tablewe show the same 
intervals for the 10 * Mq subhalo when we only use z, x, and density 
information. As we discussed in Section these observables are 
sufficient to constrain the subhalo properties up to the Id degen- 
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]—[ SDSS 
I—J Super-SDSS 
]—[ LSST 
I—I Error Free 
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Figure 9. Comparison of true and inferred properties for all setups. The error bars show the Icr scatter around the maximum likelihood value. Note that the 
error bars for each observational setup are shifted slightly along the x axis to allow comparison. On the top left panel we show the true subhalo mass versus 
the inferred mass. On the top right panel we show the true subhalo scale radius versus the infen'ed scale radius. On the bottom left panel we show the true 
mass versus the inferred impact parameter. On the bottom right panel we show the true mass versus the time since impact. For each subhalo setup, we show 
the constraints for the four observational setups. From left to right, we show the SDSS errors, the Super-SDSS errors, the LSST errors, and the setup with no 
observational errors. Note that the error free setup shows some bias in b and t. This is because our model does not perfectly match the data as described in 
Section im 


Survey 

M = 10’Mo, r, = 250 pc 

M = 10’ ®Mo, rj = 395.3 pc 

M = 10**Mo, rx = 625 pc 

SDSS 

M 6 [4.49 X 10®Mo, 1.94 x IO’Mq] 
rj 6 [53pc, 300pc] 

Me [1.27 X 10’Mo, 4.66 X IO’Mq] 
rx e [193pc, 358pc] 

M £ [6.76 X 10’Mo, 2.05 x 10**Mo] 

Tj £ [495pc, 687pc] 

Super-SDSS 

M £ [4.69 X 10®Mo, 1.85 x IO’Mo] 
f’s 6 [171pc, 341pc] 

Me [1.46 X 10’Mo, 5.22 X IO’Mq] 
rx e [380pc,435pc] 

M £ [4.30 X 10’Mo, 1.49 x 10**Mo] 

Tj £ [580pc, 634pc] 

LSST 

M 6 [5.09 X 10®Mo, 1.60 x IO’Mo] 
ts 6 [217pc, 296pc] 

M £ [1.39 X 10’Mo, 4.73 x IO’Mq] 
rx e [353pc, 397pc] 

M £ [4.27 X 10’Mo, 1.56 x 10**Mo] 

Tj £ [603pc, 631pc] 


Table 2. MCMC posteriors for vaiious survey configurations. For each survey and subhalo we list the 5-95% confidence interval of the posterior. Note that 
the mass en'ors do not change appreciably as we change the survey type while the scale radius eiTors show a more dramatic difference. This is due to the 
degeneracy between mass and velocity as we discussed in Sec.|6.4| 


eracy between mass and velocity. As expected, the constraints are 
similar to when we use all of the information. 

Figure 0 presents the posterior distribution for the IO’Mq 
subhalo when the LSST errors are used. This figure can be com¬ 
pared against Figure where we showed the posteriors with no 
observational error. While there is more scatter for the LSST er¬ 
rors, the overall trend is similar. Most of the parameters are well 
constrained except for the degeneracy between M and w. There is 


a hint of the degeneracy between b and Tj where b~ + r] is constant 
and Wx = 0 which we mentioned in Section]^ This can best be seen 
in the 2d posteriors of b - Wx and b - r,. This degeneracy will be 
much more vivid when we use SDSS errors. 

Figure [T^ demonstrates the posterior distribution for the 
IO’Mq subhalo when the SDSS errors are used. This setup has sig¬ 
nificantly more scatter than the LSST setup. The degeneracy be¬ 
tween M and w is present as expected. In addition, the degeneracy 
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Figure 10. Comparison of tme and infen'ed properties for the lO^M© subhalo with varying amounts of observables using LSST errors. The error bars show the 
Icr scatter around the maximum likelihood value. Note that the error bars for each obseiwational setup are shifted slightly along the x axis to allow comparison. 
On the left panel we show the inferred subhalo mass and on the right panel we show the inferred scale radius. We show four setups with varying amounts of 
information. First, on the left (in blue), we have “All obs” which uses all of the observables. Moving to the right (in red) we next have “No p” which uses all 
of the observables except for density. Next on the right (in green) is the setup which uses z, x, and p. Finally, the rightmost setup (in black) uses only z and p. 
We see that the constraints are similar for the first three setups. As we discussed in Section]^ these three setups have enough constraining power to constrain 
all of the observables. If we only measure z and p an additional degeneracy is possible and we see a larger error bar in r^. 




i’C ) 


Figure 7. Observables for lO^M© subhalo with Super-SDSS errors. The red 
curves show the maximum likelihood result and the blue error bars show the 
Icr eiTors of the mock observation. 


Survey 

M = 10*Mo, r, = 625 pc 

SDSS 

M E [6.76 X 10’Mo, 2.05 x 10*Mo] 
rj 6 [495pc, 687pc] 

Super-SDSS 

M E [4.37 X 10’Mo, 1.54 x 10*Mo] 

Ts 6 [571pc, 642pc] 

LSST 

M E [4.56 X 10’Mo, 1.56 x 10*Mo] 

Ts 6 [608pc, 640pc] 


Table 3. MCMC posteriors for various survey configurations when only 
z,x, and p are measured. For each survey and subhalo we list the 5-95% 
confidence interval of the posterior. Since these three observables are as 
constraining as when we use all observables (see Sec. |5.3.1) , these ranges 
are very similar to those reported in Table[^ 


Figure 8. Observables for 10^M© subhalo with SDSS errors. The red curves 
show the maximum likelihood result and the blue en'or bars show the Icr 
errors of the mock observation. Note that there are no error bai's for Ax, Vy, 
and x’z since these are not used for fitting. 


between and b when = 0 is now evident. It is especially clear 
in the 2d posterior of b - r,,. This degeneracy is due to the fact that 
the SDSS setup has no leverage on the distance to the stream and 
only a modest constraint on the radial velocity. As we discussed in 
Section]^ this allows for an additional degeneracy to appear. Note 
that this is not a true degeneracy since we have some constraint on 
the radial velocity and the gap size has a weak dependence on b. 
As a result, the likelihood is still peaked near the true value of 
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M (10^ Mq) (km/s) vjy (km/s) w, (km/s) r, (kpc) b (kpc) t (Myr) 


Figure 11. Result of MCMC for lO^M© subhalo with = 250 pc with LSST errors. In addition to the degeneracy between M and w, there is a hint of the 
degeneracy between b and which we describe in the text. This degeneracy is much more apparent in Fig. where we use SDSS errors. Apart from these 
degeneracies, most of the parameters ai'e well constrained in this setup. 


7 MASS RANGE PROBED BY CUIfRENT AND FUTURE 
SURVEYS 


In Section]^ we showed the results of a suite of MCMCs where 
we were able to measure the properties of subhaloes down to Af ~ 
IO’Mq. While this explicit demonstration is useful, we can also use 
the analytic results of Section to build some intuition as to the 
constraining power of each observable. We will restrict the analysis 
to the case of a direct impact, i.e. b = 0, since these create the 
largest signature. The detectability of a gap is controlled by the 
maximum kick imparted by the subhalo. This maximum kick is 
given by the maxima of Q: 


max|Av;r| = 
maxIAVj,] = -^ 

max I A Vj. I = 


GM wiiWx 
wrs ww± ’ 
GM 

wrs w 

GM W||Wj 
wr, ww± 


Thus we see that the maximum velocity kicks are roughly 
. GM 

Av Si -. 


(25) 


(26) 


We can now plug this expression into those for each of our observ¬ 
ables and compare the maximum deviations of the stream against a 
typical uncertainty on the centroid from the LSST setup. We note 


© 0000 RAS, MNRAS 000, 000-000 








































Properties of Dark Subhaloes from Gaps in Tidal Streams 13 



E 


E 


E 


u 

Q. 


U 

Q. 































o V n.- 'b- ^ 


0 > 




M (10^ M„) w^ (km/s) 


Wj, (km/s) 


w^ (km/s) 


r, (kpc) 


b (kpc) 


t (Myr) 


Figure 12. Result of MCMC for IO’Mq subhalo with Tj = 250 pc with SDSS errors. In addition to the degeneracy between M and w, there is also a degeneracy 
between b and rj where + d is constant. This degeneracy occurs for Wx = 0 and is discussed in the text. 


that these errors are significantly less lhan the measurement error 
for a single star. 

7.1 xysij/ 

The maximum radial deflection is given by the maximum of § 
which gives 


where w is the relative velocity of the flyby so we see that lower 
velocity flybys allow us to probe a lower mass range. We have re¬ 
arranged this expression to isolate the parameters describing the 
subhalo since these are the most relevant for understanding Galactic 
substructure. For typical values of w, v^, ~ 200 km/s, ro = 10 kpc, 
y- = 2 (i.e. a flat rotation curve), and a radial distance error of 
Axenor = 15 pc, we get 


Ay — 
^-A^max 


2roAv 


(i + 


(27) 


This must be larger than the error in estimating the stream centroid, 
AjCcrror. To determine the mass range accessible with a given error, 
we plug in the expression for Av, i.e. j26|>, and get 


M w V,y^AXe,.„r 

^ 2ro{l + 


(28) 


W ^0 

200 km/s / \ lOkpc 

(29) 

This expression can be used to understand the mass range probed as 
we change the physical and observational setup. At the edge of this 
limit, i.e. if M/r, = 6 x 10*Mo/kpc, we should expect a signal-to- 
noise ratio of 1. Note that we have not included the y dependence 
in this expression. 


M Ax.™, 

— >6x 10'’-^ —^ 
Tj kpc \ 15pc 


200 km/s 
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7.2 zysif/ 


The maximum deflection perpendicular to the plane is given by 

Av 

A(5,,„ax = —. (30) 

This must be larger than the error in estimating the stream’s posi¬ 
tion on the sky, A^.,error- This gives 

M w 

— > 7;Vj,A(5.,eiTOr- (31) 

tr 

For the typical values of w, Vy ~ 200 km/s, with the error in the 
stream latitude is given by the error when estimating the centroid, 
Adj,error ^ 0.02°, we get 


M 

Ti 


>4x10'’ 


kpc 


/ Ad, error 
\ 0 . 02 ° 


w W \ 

200 km/s/\200km/s/' 


(32) 


For typical values of w ~ 200 km/s, = 2, and Azerror = 2 km/s at 
10 kpc, this gives 


M 


>9x10^^ 

kpc 


/ AZerror \ 

\2km/s / 


w \ 
200 km/sj 


(41) 


By comparing these results, we see that the most discerning 
observable is the shape of the stream on the sky, followed by the 
distance to the stream, the radial velocity, and finally the proper 
motions. The proper motion perpendicular to the orbital plane is 
less constraining since the velocity oscillations are larger along the 
orbit. It is not possible to make such a simple estimate for the con¬ 
straining power of the density along the stream. Note that these 
formulae can be used to quickly evaluate the mass and size range 
probed by any survey. 


7.3 xysif/ 

The maximum deflection in radial velocity is given by 

2Av I ^ 

Aw = —yi + ^. (33) 

This must be larger than ATerro,, which gives 

M ^ W yAXerror .... 

^ /-• I -■ 

For typical values with w ~ 200 km/s, y^ = 2, and ATerro, = 1 km/s, 
this gives 


M 

rs 


> 3 X 10' 


kpc 


/ ATerror \ 

\ Ikm/s / 


w \ 
200 km/s / 


(35) 


7.4 y\sif/ 

The maximum deflection in proper motion along the streams’s orbit 
is given by 



This must be larger than Ayei-ror, giving 



(36) 


(37) 


For typical values with w ~ 200 km/s, y^ = 2, and Ayenor = 2 km/s 
at 10 kpc, this gives 


y >3xl0^^f^^lf--]. 

kpc \2km/s/ \ 200km/s/ 


(38) 


7.5 z vs i/' 

The maximum deflection in the proper motion perpendicular to the 
orbital plane is given by 

Azmiix = Av. (39) 

This must be larger than Azenor, which gives 
M w 

— > (40) 

r, G 


8 DISCUSSION 

The results above show that with current and upcoming surveys, 
it is possible to infer the properties of subhaloes down to a mass 
of 10 ^ Mq at a distance of 10 kpc (note, however, that according to 
Section [53| this measurement can be extended to much larger dis¬ 
tances). In addition, we quantify the accuracy with which the ge¬ 
ometry and the dynamics of the encounter as well as the mass and 
the size of the perturber can be measured. Importantly, we demon¬ 
strate that while there exists a degeneracy between subhalo’s mass 
and velocity, its detrimental effect can be remedied with the use of 
reasonable priors. We also presented expressions for the mass range 
accessible for a given observational error in Section|^ 

These mass limits should be considered as a best case sce¬ 
nario. In reality, there will be several additional complications 
which we will now discuss. These will introduce additional errors 
which may raise the mass limit measurable with the method dis¬ 
cussed here. We will also discuss other causes of substructure in 
tidal streams and how these can be distinguished from gaps caused 
by subhaloes. 


8.1 Discrete Degeneracies 

In Section|^ we characterized the degeneracies when inferring sub¬ 
halo properties. We showed that if we use all of the observables, 
there is a continuous degeneracy between M and w where M/w is 
fixed. In addition to this continuous degeneracy, it is possible that 
there are discrete degeneracies with different points in the parame¬ 
ter space giving rise to the same gap shape. Since the gap shape os¬ 
cillates in time, if we wait roughly an orbital period, the gap shape 
will look similar in z vs ifr. Of course, the gap will now be larger so 
we must change the other parameters to get the same gap size and 
shape. Thus, we might imagine a discrete degeneracy is possible. 

We have searched for this discrete degeneracy by running 
gradient descents in our negative likelihood space and initiating 
MCMCs near the location of different likelihood peaks. In the case 
of no observational errors, we have found no evidence for such de¬ 
generacy since the log-likelihoods of the varying peaks differ by 
more than 1000 from the true peak. As a result, we do not be¬ 
lieve this is a true degeneracy. However, if we include observa¬ 
tional errors, discrete degeneracies are possible in some cases. For 
the 10 *Mq subhalo and using SDSS observations, we cannot distin¬ 
guish between solutions which differ by an orbital period. However, 
if we use Super-SDSS or LSST observations, these solutions dif¬ 
fer by roughly 5 and 10 respectively in log-likelihood. For the less 
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massive subhaloes, i.e. 10 ’-^Mq and 10’Mg, the discrete degenera¬ 
cies cannot be distinguished even with LSST quality observations. 
In this work we have only considered one setup: a direct impact on 
the stream which we observe when the oscillation in the shape on 
the sky is large. It is not clear how these effective degeneracies will 
change for a different impact geometry, if we observed the gap at a 
different time, or if the progenitor is on an eccentric orhit. We will 
explore this issue in future work. 


8.2 Eccentric orbits 

The analysis done in this work was for progenitors on circular or¬ 
bits. In reality, progenitors will be on orbits with a range of eccen¬ 
tricities. This will result in the gap stretching and compressing as 
it passes near pericenter and apocenter respectively. Since analytic 
orbits are not known for an arbitrary spherical potential, it is not 
possible to repeat the analysis presented here for eccentric orbits. 
However, the procedure outlined in this work can be done numeri¬ 
cally and the qualitative picture presented here should be the same. 
In future work, we plan to address this by re-casting this analysis 
in action-angle formalism. 


8.3 Uncertainty in the potential 

In the method above, we have assumed perfect knowledge of the 
potential. Any uncertainty in the potential, and thus the shape of 
the stream, will complicate the inference of the subhalo properties. 
However, since the potential is expected to be smooth, its effect on 
the uncertainty of properties derived from kinks and wiggles over 
~ 10° scales in the stream should not be very large. These wiggles 
can only be created by small scale structure. 


8.4 Oscillations 

One additional complication which we have mentioned above is the 
fact that all of these observables oscillate in time. Thus, the quality 
of the inference depends on what phase of the oscillation the gap 
is in. In the worst case scenario, the observables which have the 
most constraining powers (i.e. the shape on the sky or the radial 
velocity) will give less information than expected since they are in 
a phase with little signal. This will lead to a broader posterior and 
may change the results of Section so that, in a particular phase, 
the proper motion observables are more informative than the shape 
of the gap on the sky. 


8.5 NEW Subhaloes 

In this analysis, we have assumed that subhaloes can be modelled as 
Plummer spheres. In reality, the subhaloes will be closer to NFW 
profiles ( [Navarro, Frenk & White|[l997) . Since the kick induced 
by the passage of an NFW subhalo cannot be computed analyti¬ 
cally, this means we will have to compute the kicks numerically. 
However, since numerics will already be necessary because of the 
non-circular orbits, this should not introduce much additional com¬ 
plication. Furthermore, since the NFW profile is controlled by two 
parameters (concentration and virial mass), this more realistic setup 
will have the same dimensionality as in this work and should ex¬ 
hibit the same qualitative features seen in this work. 


8.6 Epicyclic Feathering 

One physical process can produce density variations in a stream 
is epicyclic feathering (e.g. Kiipper, MacLeod & Heggie| |2008] 
[Ktipper et al.||20i0) l. In order to carry out the process described 
in this work on real streams, we will need to start with a detailed 
stream model. This model will include these epicyclic feathers so 
we will be able to determine the likelihood of a gap feature having 
arisen from the epicyclic feathering. More generally, the epicyclic 
feathers have different features from the gaps: the epicycles resem¬ 
ble cycloids while the stream near a gap looks like a sinusoid. Thus, 
with sufficient observations, it should be possible to distinguish the 
two. This will be easiest if the gap is far from the progenitor, prob¬ 
ing an older part of the stream where it is well-ordered. Fortunately, 
these older parts of the stream are also the most likely to have gaps 
since they have been around the longest. On a promising note, the 
N-body simulations in our work naturally included this epicyclic 
feathering and our model was robust enough to make accurate pre¬ 
dictions despite not accounting for these effects. 


8.7 Gap structure from other substructure 

In addition to dark substructure in the Galaxy, there are also bary- 
onic structures like giant molecular clouds and spiral waves. These 
can potentially cause confusion but as we have shown in Section]^ 
and Section]^ the time since impact can be well constrained. This 
allows us to rewind the position of the gap back to where the impact 
occurred. If this results in the gap being created in or very close to 
the disk, baryonic effects should be considered. However, if the im¬ 
pact occurs far from the disk, this should allow us to rule out such 
effects. In addition, the results of this work show it is possible to get 
independent constraints on the mass and scale radius. This should 
give even more evidence to rule out baryonic structures if the gap 
is created from a large and low-density structure. 


8.8 Searching for Stream Gaps 

The primary signatures of gaps are the oscillations in the stream on 
the sky outlined in this work and the contrast between the low cen¬ 
tral density and enhanced density on the edge of the gap discussed 
in Paper I. The results of this work and the results of Paper I can be 
used to create match filters to search for these features. However, 
a simpler approach may be to search for qualitative features in the 
data. Namely, an underdensity with overdensities around it, and an 
oscillation in the shape of the stream on the sky. This requires only 
positions of stars from photometric survey data and does not re¬ 
quire any velocities. Once such a feature is found, follow-up would 
be needed to constrain the radial velocities. As we discussed in 
Section |5.3.1| these three measurements would be sufficient to con¬ 
strain the subhalo properties. 


9 CONCLUSIONS 

In this work, we have shown for the first time, how the properties 
of a dark subhalo can be extracted from the measurements of the 
gap it creates in a stellar tidal stream. We used the analytic results 
of Paper I to forecast the evolution of the 6d structure of the gap 
shape arising from a subhalo flyby. We then used these predictions 
to model gaps produced in N-body simulations. We carried out in¬ 
ference in an error-free case as well as for observables with realistic 
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noise for a range of plausible instrumental setups. The main con¬ 
clusions of our study are as follows: 

(i) The information content of the subhalo-stream interaction 
is encoded in the velocity kicks received by the stream stars. 

(ii) Because the velocity kick is an integral of the acceleration 
imparted by the subhalo onto stream particles, it depends both on 
the strength of the perturber’s pull as well as the duration of the 
encounter. In other words, there exists a degeneracy between the 
mass of the subhalo and the relative velocity of the flyby. 

(iii) Using the shape of the gap without observational errors, 
it is possible to uniquely infer the impact geometry, its dynamics as 
well as the age for subhaloes with masses > 10 ^Mq. 

(iv) Importantly, it is also possible to independently infer the 
properties of the dark perturber, such as its scale length and its 
mass. 

(v) If we account for observational errors from ongoing and 
planned surveys, it is possible to infer the impact geometry, its dy¬ 
namics as well as the age for subhaloes with masses > 10 ^Mq. For 
example, for the IO’Mq subhalo with a scale radius of 250 pc, the 
Icr uncertainty on the measurement of the scale radius, the impact 
parameter, the time since flyby decreases from 73pc, 135pc, 34.8 
Myr for SDSS to 24 pc, 56 pc and 10.7 Myr for LSST. 

(vi) Due to the mass-velocity degeneracy mentioned above, 
while it is possible to infer the mass of the perturber, the accuracy 
of the measurement does not scale simply with the signal-to-noise 
ratio of the observables. Yet, in our experiments we found the un¬ 
certainty of the perturber mass to be typically on the order of 0.15 
dex. 

(vii) Since many of the gap properties oscillate in time, a very 
similar gap will be produced by subhalo hits separated by an orbital 
period. While these different solutions can be ruled out in the error- 
free case, if we add observational errors these multiple solutions 
cannot be distinguished for the IO’Mq and IO^ ^Mq subhaloes. 

(viii) Using the analytic model for the stream, we explore the 
information contained in several dilferent combinations of observ¬ 
ables. We find that just three measurements, the stream density, the 
stream shape on the sky, and the radial velocity along the stream 
contain enough information to constrain the subhalo properties. We 
believe, therefore, that the stream gap analysis presented here can 
be carried out at heliocentric distances much larger than presented 
here. 

The analytic nature of the results in this work allows us to 
build a simple intuitive picture of how to infer subhalo proper¬ 
ties from stream gaps. We discussed the difficulties which could 
arise when extending this model to gaps created in more realistic 
streams. This will require numerical modelling but we expect the 
same qualitative picture we have described here. It should be pos¬ 
sible to constrain subhaloes which impact streams on non-circular 
orbits up to the same Id degeneracy seen here. In future work, we 
will explore this by re-casting this framework in angle-action coor¬ 
dinates. 

Given the promising results of this study, the outlook for dark 
subhalo detection in the Milky Way is bright. If Dark Matter exists, 
we expect that these elusive clumps will soon be unveiled. 
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APPENDIX A: COEFFICIENTS 

For completeness, we list the coefficients Aj(p) and S,(^ which we 
introduced in Section[5T]and Section[Y2|respectively: 
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